METHOD FOR BACKUP DUAL-FREQUENCY NAVIGATION DURING 
BRIEF PERIODS WHEN MEASUREMENT DATA IS UNAVAILABLE ON 

ONE OF TWO FREQUENCIES 

[0001] The present invention relates generally to technologies associated with 

positioning and navigation using satellites, and more particularly to dual-frequency 
navigation using the global positioning system (GPS). 

BACKGROUND 

[0002] The global positioning system (GPS) uses satellites in space to locate objects 

on earth. With GPS, signals from the satellites arrive at a GPS receiver and are used to 
determine the position of the GPS receiver. Currently, two types of GPS measurements 
corresponding to each correlator channel with a locked GPS satellite signal are available for 
civilian GPS receivers. The two types of GPS measurements are pseudorange, and integrated 
carrier phase for two carrier signals, LI and L2, with frequencies of 1.5754 GHz and 1.2276 
GHz, or wavelengths of 0.1903 m and 0.2442 m, respectively. The pseudorange measurement 
(or code measurement) is a basic GPS observable that all types of GPS receivers can make. It 
utilizes the C/A or P codes modulated onto the carrier signals. The measurement records the 
apparent time taken for the relevant code to travel from the satellite to the receiver, i.e., the 
time the signal arrives at the receiver according to the receiver clock minus the time the 
signal left the satellite according to the satellite clock. 

[0003] The carrier phase measurement is obtained by integrating a 

reconstructed carrier of the signal as it arrives at the receiver. Thus, the carrier phase 
measurement is also a measure of a transit time difference as determined by the time the 
signal left the satellite according to the satellite clock and the time it arrives at the receiver 
according to the receiver clock. However, because an initial number of whole cycles in transit 
between the satellite and the receiver when the receiver starts tracking the carrier phase of 
the signal is usually not known, the transit time difference may be in error by multiple carrier 
cycles, i.e., there is a whole-cycle ambiguity in the carrier phase measurement. 
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[0004] With the GPS measurements available, the range or distance between a 

GPS receiver and each of a multitude of satellites is calculated by multiplying a signal's 
travel time by the speed of light. These ranges are usually referred to as pseudoranges (false 
ranges) because the receiver clock generally has a significant time error, which causes a 
common bias in the measured range. In addition, several error factors exist that can lead to 
errors or noise in the calculated range, such as the ephemeris error, satellite clock timing 
error, atmospheric effects, receiver npise and multipath error. The common bias from 
receiver clock error is usually solved for along with the position coordinates of the receiver as 
part of the normal navigation computation. 

[0005] With standalone GPS navigation, where a user with a GPS receiver 

obtains code and/or carrier-phase ranges with respect to a plurality of satellites in view, 
without consulting with any reference station, the user is very limited in ways to reduce the 
errors or noises in the ranges. To eliminate or reduce some of these errors, differential 
techniques are typically used in GPS applications. Differential GPS (DGPS) operations 
typically involve one or more reference GPS receivers in fixed locations, a user (or 
navigation) GPS receiver, and communication links among the user and reference receivers. 
The reference receivers are used to generate corrections associated with some or all of the 
above error factors. The corrections are supplied to the user receiver and the user receiver 
then uses the corrections to appropriately correct its computed position. 

[0006] A number of different techniques have been developed to obtain high- 

accuracy differential navigation using the GPS carrier-phase measurements. The highest 
accuracy technique is generally referred to as "real-time kinematic" (RTK) and has a typical 
accuracy of about one-centimeter. However, in order to obtain that accuracy, the whole-cycle 
ambiguity in the differential carrier-phase measurements must be determined. When the 
reference receiver is a substantial distance (more than a few tens of kilometers) from the 
navigation receiver it may become impossible to determine the whole-cycle ambiguity and 
the normal RTK accuracy cannot be achieved. Under these adverse circumstances the best 
that can be done is often to estimate the whole-cycle ambiguities as a real-valued (non- 
integer) variable. This practice is often referred to as determining a "floating ambiguity" 
value. 
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[0007] One method for determining the "floating ambiguity" value is to form 

refraction corrected code and carrier-phase measurements, scale the refraction corrected 
carrier-phase measurement to the same unit as the refraction corrected code measurement, 
and form an offset by subtracting the refraction corrected carrier-phase measurement from the 
refraction-corrected code measurement. This offset value can be recursively averaged over 
time so that it becomes an increasingly accurate estimate of the "floating ambiguity." Exactly 
the same net result can be obtained by smoothing a code measurement with a linear 
combination of the corresponding LI and L2 carrier-phase measurements that is formed to 
match the ionospheric refraction effect of the code measurement. 

[0008] Several types of differential GPS systems that provide measurements 

or measurement corrections to navigation receivers are currently available. Among them, the 
High Accuracy Nationwide Differential GPS System (HA-ND GPS), which is developed 
cooperatively by several U.S. government organizations, uses ground based reference sites. 
This system transmits the corrections to the user using Coast Guard beacons that can reach 
users at ranges of a few hundred kilometers. John Deere has developed the StarFire™ system, 
which transmits corrections via communication satellites with both a regional wide area 
correction data stream and a global DGPS correction data stream. In these systems, 
navigation results in the 10 centimeter range can be obtained after the carrier-phase floating 
ambiguities have been determined with sufficient accuracy, that is, after sufficient time has 
elapsed since the navigation receiver starts tracking the satellite signals. 

[0009] One of the principal problems of these navigation systems is that 

anything such as interfering signals, shading or signal blockage, etc., which causes one of the 
signals from any of the satellites to be temporarily lost, will cause "cycle slips" in the carrier- 
phase measurements and the floating ambiguity value will no longer be correct. In the current 
commercial environment, the L2 signals are much more apt to be lost than the LI 
measurements. There are several reasons for this. First the broadcast LI signal is stronger 
than the broadcast L2 signal. In addition, commercial access to the L2 signal requires a 
"codeless" or "semi-codeless" technique to be employed to avoid the selective availability 
imposed on the L2 signal by the military. As a result, only a small amount of interference or 
signal blockage can cause a loss of the L2 measurements. Without some means of 
reinitializing the floating ambiguity value, a long time interval will be required to determine 
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anew the correct floating ambiguity value after the L2 signal returns. Therefore there is a 
need for a technique to reinitialize the floating ambiguity value after a brief L2 signal outage 
so that the long initialization process can be avoided. 

SUMMARY 

[0010] The present invention includes a method for performing backup dual- 

frequency navigation whereby the L2 code and carrier-phase measurements are synthesized 
using a combination of the retained LI carrier-phase measurements and a model of the 
ionospheric refraction effects, which is updated when measurements on both the LI and L2 
frequencies are available. As an optional process, a divergence between the retained code and 
carrier phase measurements can be used to detect slowly changing deviations from the 
ionospheric refraction model. This allows an increase in the interval over which synthesized 
measurements can be successfully generated. 

[0011] In one embodiment of the present invention, the backup dual-frequency 

navigation is performed for each satellite from which the L2 measurements are lost for a time 
period at the user GPS receiver, and the method for performing the backup dual-frequency 
navigation includes steady-state processing when measurements on both the LI and L2 
frequencies from the satellite are available. During the steady-state processing, smoothed 
code measurements and smoothed offsets between code and carrier-phase measurements are 
computed. Also, corrections to an ionospheric model are generated. Thereafter, when direct 
measurements on the L2 frequency from the satellite are unavailable, backup operations are 
performed for each measurement epoch until the L2 signals are detected again at the user 
GPS receiver. During the backup operations, the ionospheric model corrections are used to 
generate estimated L2 carrier-phase measurements, which are used to generate estimated 
code measurements on both the LI and the L2 frequencies. The estimated and measured code 
measurements on the LI frequency are used in an optional step in which ionospheric model 
corrections are updated. Upon the return of the L2 signals, a transition to dual frequency 
navigation using both the LI and L2 signals from the satellite is performed. 
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[0012] Thus, the method in one embodiment of the present invention allows dual 

frequency operation at a GPS receiver to continue in the situation when signals from one or 
more satellites on one of the frequencies become unavailable for a time period. 

DRAWINGS 

[0013] FIG. 1 is a block diagram of a computer system that can be used to perform 

the backup dual frequency navigation method according to one embodiment of the present 
invention. 

[0014] FIG. 2 is a flowchart illustrating the method for backup dual frequency 

navigation according to one embodiment of the present invention. 

[0015] FIG. 3 is a flowchart illustrating a step for generating smoothed code 

measurements and smoothed offsets between the code and carrier-phase measurements 
during steady state processing in the method for backup dual -frequency navigation. 

[0016] FIG. 4 is a flowchart illustrating a step for generating ionospheric model 

corrections during steady state processing in the method for backup dual frequency 
navigation. 

[0017] FIG. 5 is a flowchart illustrating a step for generating synthesized (or 

estimated) L2 carrier-phase measurement in the method for backup dual-frequency 
navigation when direct L2 measurements are unavailable. 

[0018] FIG. 6 is a flowchart illustrating a step for generating synthesized code 

measurement in the method for backup dual-frequency navigation when L2 measurements 
are unavailable. 

[0019] FIG. 7 is a flowchart illustrating an optional step for updating the ionospheric 

model corrections in the method for backup dual frequency navigation when L2 
measurements are unavailable. 

[0020] FIG. 8 is a flowchart illustrating a transition to steady-state dual-frequency 

navigation after the L2 signal returns. 
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DESCRIPTION 

[0021] FIG. 1 illustrates a system 100 for performing backup dual-frequency 

navigation in case of an occasional loss-of-lock on the L2 signal from one of the satellites, 
according to one embodiment of the present invention. As shown in FIG. 1, system 100 can 
be a microprocessor-based computer system 100 coupled to a GPS receiver 110, which 
provides raw GPS observables to system 100 for processing. These observables include GPS 
code and carrier phase measurements, ephemerides, and other information obtained according 
to signals received from a plurality of satellites 101. 

[0022] To facilitate differential operations, system 100 may also be coupled to a 

reference station 120 via a radio link 124. The reference station 120 provides GPS 
observables measured thereat and/or GPS corrections calculated thereat. In wide-area or 
global applications, system 100 may be coupled to one or more central hubs 130 in 
communication with a group of reference stations (not shown) via radio and/or satellite links 
134. The hub(s) 130 receives GPS observables from the group of reference stations and 
computes corrections that are communicated to the system 100. 

[0023] In one embodiment of the present invention, system 100 includes a central 

processing unit (CPU) 140, a memory device 148, a plurality of input ports 153, 154, and 
155, one or more output ports 156, and an optional user interface 158, interconnected by one 
or more communication buses 152. Memory 148 may include high-speed random access 
memory and may include nonvolatile mass storage, such as one or more magnetic disk 
storage devices. Memory 148 may also include mass storage that is remotely located from 
the central processing unit 140. Memory 148 preferably stores an operating system 162, a 
database 170, and GPS application programs or procedures 164, including procedures for 
backup dual frequency navigation 166 according to one embodiment of the present invention. 
The operating system 162 and application programs and procedures 164 stored in memory 
148 are for execution by the CPU 140 of the computer system 100. Memory 148 preferably 
also stores data structures used during the execution of the GPS application procedures 166, 
such as GPS measurements and corrections, as well as other data structures discussed in this 
document. 
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[0024] The input ports 154 are for receiving data from the GPS receiver 1 10, the 

reference station 120, and/or the hub 130, respectively, and the output port(s) 156 can be used 
for outputting calculation results. Alternately, calculation results may be shown on a display 
device of the user interface 158. 

[0025] The operating system 162 may be, but is not limited to, the embedded 

operating system, UNIX, Solaris, or Windows 95, 98, NT 4.0, 2000 or XP. More generally, 
operating system 162 has procedures and instructions for communicating, processing, 
accessing, storing and searching data. 

[0026] As indicated by the dashed line 105 in FIG. 1, in some embodiments, the GPS 

receiver 110 and part or all of the computer system 100 are integrated into a single device, 
within a single housing, such as a portable, handheld or even wearable position tracking 
device, or a vehicle-mounted or otherwise mobile positioning and/or navigation system. In 
other embodiments, the GPS receiver 1 10 and the computer system 100 are not integrated 
into a single device. 

[0027] FIG. 2 is a flowchart illustrating a process 200 for performing backup dual- 

frequency navigation according to one embodiment of the present invention. The process 200 
is performed for each satellite 101 from which the L2 measurements are lost for a time period 
at the GPS receiver 1 10. As shown in FIG. 2, process 200 includes steps 210 and 220, which 
are performed during steady-state processing when measurements on both the LI and L2 
frequencies from the satellite are available. In step 210, smoothed code measurements and 
smoothed offsets between code and carrier-phase measurements are computed. In step 220, 
ionospheric model corrections are generated. Thereafter, when direct measurement on L2 
frequency from the satellite becomes unavailable, steps 230, 240, and optional step 250 are 
performed for each measurement epoch before the L2 signals returns at the GPS receiver 1 10. 
In step 230, the ionospheric model corrections are used to generate estimated L2 carrier- 
phase measurements, which are used in the subsequent step 240 to generate estimated code 
measurements on both LI and L2 frequencies. The estimated and measured code 
measurements on the LI frequency are used in the subsequent optional step 250 in which 
ionospheric model corrections are updated. The process 200 then proceeds to a step 260 in 
which it is determined whether L2 signals from the satellite have returned. If L2 signals have 
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not returned, steps 230 through 250 are repeated for the next measurement epoch using the 
updated ionospheric model corrections. Otherwise, upon the return of L2 signals, a transition 
to dual frequency navigation using both LI and L2 signals from the satellite is performed in 
step 270. 

[0028] During steady-state processing when measurements from both LI and L2 

frequencies are available, the multipath error in each code measurement can be minimized by 
forming a combination of the LI and L2 carrier-phase measurements that matches the 
ionospheric refraction effect in the code measurement, and by smoothing the code 
measurement with the carrier-phase measurement combination. Many receivers make both a 
C/A-code measurement and a P-code measurement on the LI frequency. Either of the C/A or 
P-code measurement can be used as the LI code measurement. However, whichever of the 
two is chosen, the same should be used at the user and the reference station(s) since small 
biases exist between the two measurements. In the discussion that follows, the LI frequency 
(equal to about 1.57542 GHz) is designated as f\ and the L2 frequency (normally equal to 
about 1.2276 GHz) is designated as/ 2 . The pseudorange code measurement (whether C/A or 
P) on the LI frequency is designated as Pi and the pseudorange code measurement on the L2 
frequency is designated as P 2 . The LI carrier-phase measurement in meters will be designated 
simply as Li and the L2 carrier-phase measurement in meters will be designated as L 2 . The 
carrier-phase measurements are scaled by the wavelengths and an approximate whole-cycle 
ambiguity value is added to each so that the phase measurements are made close to the same 
value as the corresponding code measurement. Thus, using 0? to designate the raw phase 
measurement in cycles at the/i frequency and (fc to designate the raw phase measurement in 
cycles at the f 2 frequency, we have the following relationships: 

L, ={<t> x +N»)A x (1) 

L 2 =(0 2 +N° 2 )A 2 (2) 

[0029] The wavelength X } for the LI frequency is approximately equal to .1903 

meters and the wavelength of X 2 for the L2 frequency is approximately .2442 meters. The 
approximate whole-cycle values of, and N 2 are added at the start of carrier-phase 
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tracking to give values that are within one wavelength of the corresponding code 
measurements simply to keep the differences to be formed subsequently small. 

[0030] FIG. 3 is a flowchart illustrating in more detail step 210 in process 200, in 

which smoothed code measurements and smoothed offsets between the code measurements 
and corresponding carrier-phase measurements are computed during steady-state processing 
when signals on both LI and L2 frequencies are available from the satellite. When the L2 
signal is not available, the previously computed values for the smoothed PI offset (Oi), 
smoothed P2 offset (0 2 ) and the estimated AAT,^, - AN 2 A 2 (O2-O1) from the last epoch of 
steady-state processing are stored and used during backup dual frequency operation. 

[0031] As shown in FIG. 3, step 210 includes a substep 310, in which a first linear 

combination Mi of Li and L 2 are formed to match the delay due to the ionospheric refraction 
effect on code measurement Pi, and a substep 320, in which a second linear combination M 2 
of Li and L 2 are formed to match the delay due to the ionospheric refraction effect on code 
measurement P 2 . Substeps 310 and 320 are performed according to the following equations: 

M x = (K i +K 2 )L l -2K 2 L 2 (3) 



M 2 =2K l L l -(K l +K 2 )L 2 (4) 
where K } and K 2 are coefficients defined as follows: 



ft 

A 2 - fi 



= - =2.5457 (5) 



K 2 = Jl =1.5457 (6) 

[0032] Because the ionospheric effects on the code measurements Pi and P 2 have 

been matched by the respective linear combinations Mi and M 2 of the carrier-phase 
measurements, and because all clock variations and motions for either the satellite transmitter 
or the user receiver have identical effects on the code and carrier-phase measurements, Mi 
and Pi, or M 2 and P 2 , should be identical except for possible whole-cycle ambiguity errors in 
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the carrier-phase combination, Mi or M 2 , and the higher multipath noise in the code 
measurement Pi or P 2 , respectively. This allows the formation of smoothed code 
measurements which approaches the small measurement noise of the carrier-phase 
measurements but without the associated whole-cycle ambiguity. 

[0033] Thus, step 210 further includes a substep 330, in which an offset between Pi 

and Mi is computed, and a substep 350, in which the offset is processed in a low pass filter to 
form a smoothed offset Oi between Pi and Mi (referred in FIG. 3 and subsequently as the 
"smoothed Pj offset"). In parallel, step 210 also includes a substep 340, in which an offset 
between P 2 and M 2 is computed, and a substep 360, in which the offset is processed in a low 
pass filter to form a smoothed offset O2 between P 2 and M 2 (referred in FIG. 3 and 
subsequently as the "smoothed P 2 offset"). Using subscript "f ' to designate the measurements 
at a specific measurement epoch, the low pass filter in substep 350 or 360 forms the 
smoothed Pi or P 2 offset by sequentially averaging the offset according to the following 
equation: 

Ojlj = 0 A i _ x + (P AJ - M A i - 0 A i _ x )ln (7) 

where X = 1 or 2 for designating the LI or L2 frequency, and O^i represents the smoothed 
Pior P 2 offset at the i th measurement epoch. The low pass filter in substep 350 or 370 forms 
sequential averages until a maximum averaging interval is achieved and then it converts to an 
exponential smoothing filter. So, n equals to i until the maximum averaging interval is 
reached and then holds at that maximum value afterwards. It should be noted that other forms 
of low-pass filtering could be used. One alternative is to model the multipath errors in the 
code measurements as correlated noise and use a stochastic model of the multipath error in a 
Kalman filter to obtain an estimated offset between the code and carrier-phase 
measurements. 

[0034] Step 210 in the process 200 further includes substeps 370 and 380, in which 

the smoothed P| and P 2 are each formed by summing the corresponding offset with the 
corresponding carrier-phase measurement, as in the following: 

S a =O a +M a (8) 
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where Sx, A, = 1 or 2, represents the smoothed Pi or P 2 code measurements. 

[0035] It is noted that the values of the smoothed Pi and P2 offsets will approach 

specific values as the number of measurement epochs used in the smoothing process (referred 
herein also as the "averaging interval" or "smoothing count") increases. Specifically, when 
enough averaging has been performed, the following should hold, 

O x = (K x + K 2 )AN l A l - 2K 2 AN 2 A 2 (9) 

0 2 = 2K X AN X A X - (K x + K 2 )AN 2 A 2 (10) 

where the values of AN\ and AN 2 represent the errors in the initial assignment A^ 0 and N 2 of 
the integer ambiguities in the raw carrier-phase measurements <t>i and fa respectively. For 
subsequent use, step 210 further includes a substep 390 in which the difference between the 
two smoothed offsets are computed to yield an estimated t±N x \ - AAf 2 A 2 : 

0 2 - O x = AN X A X - &N 2 A 2 (11) 

[0036] FIG. 4 is a flowchart illustrating in more detail the processing for generating 

ionospheric refraction corrections in step 220 in process 200. The ionospheric refraction 
corrections generated in step 220 are to be used to synthesize L2 measurements when direct 
L2 measurements are not available. As shown in FIG. 4, step 220 includes a substep 410, in 
which an ionospheric model is used to compute a modeled ionospheric bias term, I m , and 
optionally a modeled ionospheric rate term, Delta I m . The ionospheric rate term is computed 
from sequential differences of the ionospheric bias terms obtained from the model. Any of 
several ionospheric models could be used in substep 410, including the ionospheric model in 
the Wide Area Augmentation System (WAAS), whose corrections are broadcast from the 
WAAS communication satellites, the real-time ionospheric model used by the International 
GPS Service (IGS), and the ionospheric model whose corrections are broadcast from the GPS 
satellites. Since most ionospheric models generate the ionospheric refraction bias term and 
rate term in the Pi code measurement at the fi frequency, the modeled bias term and rate term 
need to be divided by the ^coefficient to obtain the expected difference between 
ionospheric delays in the Pi and P 2 code measurements. Thus, step 200 further includes a 
substep 420, in which I m and Delta I m are divided by K 2 for subsequent use. 
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[0037] Step 220 in process 200 further includes a substep 430, in which the smoothed 

code measurements computed in step 210 according to Equations (1) through (8) are 
differenced to yield a measured ionospheric bias term, and a substep 440, in which \ JK 2 is 
subtracted from the measured ionospheric bias term to produce a correction, AI , to the 
modeled ionospheric bias term. Substeps 430 and 440 are performed according to the 
following equation: 

M = S 2 -S i -IJK 2 (12) 

[0038] To generate an optional correction to the modeled ionospheric rate term, step 

220 in process 200 further includes a substep 450, in which a difference between the L2 
carrier-phase measurements taken at two consecutive measurement epochs (Delta L 2 ) is 
subtracted from a difference between the LI carrier-phase measurements taken at the two 
consecutive measurement epochs (Delta Lj) to yield a measured ionospheric rate term. 
Substep 450 is followed by a substep 460, in which (Delta l m )/K 2 is subtracted from the 

measured ionospheric rate term to produce a correction, A / , to the ionospheric rate term. 
This ionospheric rate needs to be lightly filtered to provide some smoothing without 
excessive delay. Thus, step 220 in process 200 may further include a substep 470, in which 
the result from substep 460 is processed in a low-pass filter to produce a lightly filtered 
ionospheric rate correction. This lightly filtered value of ionospheric rate correction (filtering 
equation not shown) is used subsequently in equation (15) below. By differencing the 
measured ionospheric values from the modeled values, it should be possible to generate valid 
estimates of the ionospheric effect for longer time intervals since a major portion of the 
ionospheric dynamics is handled by the model. In equation form, steps 450 to 460 can be 
represented by: 

A/ = (L M -L U-1 )-(L 2 , -L 2A _ x )-{I mi -I m ^)/K 2 (13) 

where subscript i designates the current measurement epoch, and subscript i-1 designates the 
measurement epoch prior to the current measurement epoch. 

[0039] Steps 210 and 220 in process 200, in which values such as the smoothed code 

measurements and the corrections to the ionospheric bias term and the optional rate term are 
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generated, are performed when measurements from both frequencies are available. Given that 
a sufficient interval of smoothing has occurred in the initial processing such that the values 
generated in steps 210 and 220 have most of the code multipath noise smoothed out by 
averaging, these values can be used to generate synthesized f 2 measurements in steps 230 
through 250 when measurements on the f 2 frequency are unavailable. 

[0040] FIG. 5 illustrates a process flow in step 230, in which the L2 carrier-phase 

measurement is synthesized when direct measurements on the f 2 frequency are unavailable. 
As shown in FIG. 5, step 230 in process 200 includes an optional substep 510, in which the 
correction for the ionospheric bias term generated in the previous measurement epoch and the 
modeled ionospheric bias term generated in the current measurement epoch are summed to 
produce an estimated ionospheric bias term, Inmate • Ste P 230 further includes an optional 
substep 520, in which the correction to the ionospheric rate term generated while the L2 
measurements were available is multiplied by the time period At since the L2 measurements 
became unavailable and the product of the multiplication is added to the estimated 
ionospheric bias term to produce an updated estimate of the ionospheric bias term I^ te . 

Step 230 further includes a substep 530, in which the updated estimate of the ionospheric bias 
term is subtracted from a sum of the LI carrier-phase measurement at the present 
measurement epoch and the estimated AA^ - AN 2 A 2 to produce the synthesized L2 carrier- 
phase measurement L 2 . In equation form, substeps 510, 520, and 530 can be described 
respectively by Equations (14), (15), and (16), as in the following: 



/ 



Bias 

Estimate 



= IJK 2 +AI 



(14) 



Bias 
Update 



= 1 



Bias 
Estimate 



-A/A/ 



(15) 



L 2 = L,+(AiV 2 /l 2 -A7V 1 /l 1 )-/, 



• Bias 
Update 



(16) 



where L 2 designates the synthesized L 2 . 



[0041] FIG. 6 is a flowchart illustrating in more detail the processing in step 240, in 

which the smoothed code measurements are synthesized from the LI carrier-phase 
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measurement and the synthesized L2 carrier-phase measurement. It might seem odd that the 
raw code measurement, Pi, is not used in synthesizing the smoothed code measurement at 
either frequency. Attempting to smooth the raw code measurement with the help of the 
synthesized L2 carrier-phase measurement would cause any errors in the modeled 
ionospheric refraction to generate biases that would be filtered into the offset values 
represented by equations (9), (10) and (11). To avoid creating an ionospheric refraction bias 
in the offset values, a process which is parallel to that shown in Figure 1 is used, except that 
instead of an input of the code measurements and an output of the offsets, the offsets are 
input and the synthesized code measurements are output. 

[0042] Accordingly, as shown in FIG. 5, step 240 includes a substep 610, in which 

the measured LI measurement Li and the synthesized L2 measurement L 2 are combined to 
form a carrier-phase combination M x with an ionospheric delay that matches the ionospheric 
delay in the LI code measurement Pi, and a substep 620 in which the measured LI 
measurement Li and the synthesized L2 measurement L 2 are combined to form a carrier- 
phase combination M 2 with an ionospheric delay that would match the ionospheric delay in 
the undetected L2 code measurement. In equation form, substeps 610 and 620 can be 
expressed as: 

M x = (K x + K 2 )L l - 2K 2 L 2 (17) 
M 2 = 2K l L l -(K x +K 2 )L 2 (18) 

[0043] Step 240 in process 200 further includes a substep 630, in which the smoothed 

PI offset Oi computed in step 210 is added to M x , resulting in an estimated smoothed LI 
code measurement S x , and a substep 630 in which the smoothed P2 offset 0 2 is added to M 2 , 

resulting in an estimated smoothed L2 code measurement S 2 , as expressed by the following 
equations: 

S x =M x +O x (19) 
S 2 =M 2 +0 2 (20) 
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[0044] While the raw Pi code measurement was not used to synthesize the smoothed 

code measurements, it can be used in the optional step 250 in process 200 to correct for small 
ionospheric refraction errors, which would otherwise accumulate. FIG. 7 is a flowchart 
illustrating in more detail the processing performed in the optional step 250 in process 200. 
Because the raw Pi code measurement is noisy, it must be filtered heavily in a low-pass filter 
to avoid introducing more errors from the multipath effects than it removes from ionospheric 
refraction effects. Also, because the synthesized Pi code measurement is generated from the 
LI carrier-phase measurement, any error in the ionospheric model should affect the 
synthesized Pi code measurement in a direction opposite to the way that error affects the raw 
Pi code measurement. 

[0045] Thus, step 250 includes a substep 710, in which the difference between the 

measured and synthesized code measurements is divided by 2K 2 to produce an ionospheric 
adjustment that scales with the ionospheric bias term and the optional rate term, and a substep 
720, in which this ionospheric adjustment is smoothed in a low-pass filter to remove the 
multipath errors. Step 250 further includes an optional substep 730, in which the smoothed 
ionospheric adjustment is added to the correction to the ionospheric rate term to produce an 
updated correction to the optional ionospheric rate term, and a substep 740, in which the 
smoothed ionospheric adjustment is added to the correction to the optional ionospheric bias 
term to produce an updated correction to the ionospheric bias term. 

[0046] It is also possible that a two-state estimator, e.g. an alpha-beta or Kalman 

filter, could be used to generate the updated correction to the ionospheric rate term. See Yang 
et al., "LI Backup Navigation for Dual Frequency GPS Receiver," Proceedings of the 16 th 
International Technical Meeting of the Satellite Division of the Institute of Navigation 
GPS/GNSS Conference, Sept. 9-12, 2003, Portland Oregon, which is incorporated herein by 
reference. By using some form of the process shown in FIG. 7, it may be possible to extend 
the time period that can be covered by the synthesis procedure in process 200. 

[0047] FIG. 8 is a flowchart illustrating in more detail the processing in step 270 in 

process 200, in which a transition to dual-frequency navigation is performed upon a 
determination in step 260 that the L2 signal has returned. Two tests are needed to determine 
whether or not the "floating integer" offsets computed in step 210 can be safely adjusted to 
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avoid a reinitialization of the long smoothing process otherwise required. As shown in FIG. 
8, the first test is performed in a substep 820, in which it is determined whether or not the 
interval of time At over which the L2 signal was lost exceeds a predetermined threshold. If 
the threshold is exceeded, then no adjustment is attempted and the smoothing process is 
reinitialized in a substep 830. Otherwise, the second test is performed in substeps 840 and 
850, in which the difference between the measured and the synthesized or estimated L2 
carrier-phase measurements is divided by the L2 wavelength to see if the result is close to an 
integer, i.e.: 

(L 2 - L 2 )//l 2 = integer (21) 

If the result is not within some predetermined vicinity of an integer value, substep 830 is 
performed subsequently, in which the smoothing process is reinitialized. Otherwise, the 
result is used to adjust either the floating-ambiguity in the L2 carrier-phase measurement or 
the P2 code offset value so that the code smoothing process in step 210 can be resumed after 
this simple adjustment. 

[0048] Because in practice the LI signal is virtually never lost without a concomitant 

loss of the L2 signal, the technique described herein achieves its primary intended purpose 
when used to synthesize the L2 measurements from the LI measurements during loss of only 
the L2 measurements. The present invention, however, can be applied to synthesize any of 
the LI and L2 measurements, or measurements in some other frequency, such as the L5 
frequency (equal to about 1.17645 GHz), by using measurements from another frequency that 
is not lost, with the help of a model of the ionospheric refraction effects, which is corrected 
by measurements taken while both frequencies are available. 
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